Astronomy & Astrophysics manuscript no. 12601 


© ESO 2009 


November 5, 2009 





Transneptunian objects and Centaurs from light curves 

Duffard, R. 1 , Ortiz, J.L. 1 , Thirouin, A. 1 , Santos-Sanz, P. 1 , and N. Morales 1 

1 Instituto de Astroftsica de Andaluci'a - CSIC, Apt 3004, 18080 Granada, Spain 
received/accepted 

ABSTRACT 

Aims. We compile and analyze an extended database of light curve parameters scattered in the literature to search for correlations and 
study physical properties, including internal structure constraints. 

Methods. We analyze a vast light curve database by obtaining mean rotational properties of the entire sample, determining the spin 
frequency distribution and comparing those data with a simple model based on hydrostatic equilibrium. 

Results. For the rotation periods, the mean value obtained is 6.95 h for the whole sample, 6.88 h for the Trans-neptunian objects 
(TNOs) alone and 6.75 h for the Centaurs. From Maxwellian fits to the rotational frequencies distribution the mean rotation rates are 
7.35 h for the entire sample, 7.71 h for the TNOs alone and 8.95 h for the Centaurs. These results are obtained by taking into account 
the criteria of considering a single-peak light curve for objects with amplitudes lower than 0. 15 mag and a double-peak light curve for 
objects with variability >0. 15mag. We investigate the effect of using different values other than 0.15mag for the transition threshold 
from albedo-caused light curves to shape-caused light curves. The best Maxwellian fits were obtained with the threshold between 
0.10 and 0.15 mag. The mean light-curve amplitude for the entire sample is 0.26 mag, 0.25 mag for TNOs only, and 0.26 mag for the 
Centaurs. The Period versus B-V color shows a correlation that suggests that objects with shorter rotation periods may have suffered 
more collisions than objects with larger ones. The amplitude versus H,, correlation clearly indicates that the smaller (and collisionally 
evolved) objects are more elongated than the bigger ones. 

Conclusions. From the model results, it appears that hydrostatic equilibrium can explain the statistical results of almost the entire 
sample, which means hydrostatic equilibrium is probably reached by almost all TNOs in the H range [-1,7]. This implies that for 
plausible albedos of 0.04 to 0.20, objects with diameters from 300 km to even 100 km would likely be in equilibrium. Thus, the great 
majority of objects would qualify as being dwarf planets because they would meet the hydrostatic equilibrium condition. The best 
model density corresponds to 1 100 kg/m 3 . 
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1. Introduction 

The belt of remnant planetesimals in the so-called Kuiper belt, 
whose existence was confirmed observationally in 1992 (Jewitt 
|& Luu T993J, is an important source of knowledge about the 
formation and early evolution of the Solar System. Features in 
the orbital distribution of the belt can provide important infor- 
mation about the early dynamical processes such as planetary 
migration, which was merely untested theories around 25 years 
ago (Fernandez & Ip 1984) that now can be regarded as being 
well proven. The detailed architecture of the Kuiper belt might 
also be compatible with a very unstable phase when Jupiter and 



Saturn entered into reson ance (|Tsiganis et al.| 2005 
let al _ 



Morbidelli 



J20051 |Gomes et aL|2005| ) around 800Myr after the Solar 
System formation. Although this particular dynamical scenario 
is still disputed, there is no question that the architecture of the 
transneptunian belt can offer invaluable clues about the dynam- 
ical mechanisms acting in the early Solar System. In addition 
to dynamical inferences from the orbital characteristics of the 
TNOs, much can be learnt about the physical processes that took 
place in the solar nebula by means of the study of the physi- 
cal properties of TNOs and related groups of objects. Since the 
TNOs are understood to be mixtures of ices and rocks from 
which the Centaurs and ultimately the Jupiter family comets 



(JFCs) originate (|Jewitt 2004), ground and space-based observa- 



tions of all these populations can give us clues about the physics 



of the outer Solar System. From the "Nice" dynamical model 
( fTsiganis et al.|2005[|Morbidelli et al.|2005 ~ 



Gomes et al. 2005 ), 



it is also likely that the Jupiter trojans are closely related to the 
TNOs ( |Morbidelli et al.|[2005) , but this is not yet well estab- 
lished. 

These objects are regarded as the least evolved objects in the 
Solar System, but they have experienced many different types 
of physical processes that have altered them considerably. These 
processes range from collisions, to space weathering. The col- 
lisional processes have probably left clear imprints in both the 
spin distribution and in the light-curve amplitude distribution of 
TNOs and Centaurs. The currently accepted theory based on the 
study of Davis & Farinella ( 1997 ) on collisional evolution of the 
TNOs is that objects of sizes larger than 100 km are very slightly 
collisionally evolved, whereas the smaller fragments have expe- 
rienced more collisional events and hence should have their ro- 
tations and shapes more significantly altered. This picture is also 
consistent with newer collisional evolution models (Ben avidez| 
& Campo Bagatin 2009). However, detailed modelling of the 
evolution of the spin rates related to collisional models has not 
yet been carried out for the transneptunian belt. 

Due to the dedicated efforts of two main groups studying 
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and occasional contributions by other authors, a database of 
about 40 rotation periods and 80 light-curve amplitudes could be 
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extracted from the literature in 2008. Thus, we now have suffi- 
cient numbers of TNOs light curves to start analyzing them from 
a statistical point of view. In a review by Sheppa fd et al.| p008 ), 
those data were compiled and some conclusions were outlined. 
After that compilation, a new set of light curves from 29 TNOs 
and Centaurs became available (Thiro uTn et al.|20 09) and those 
authors also presented evidence of possible biases that can be 
present in the light curve database gathered from the literature. 

We also developed a simple model to interpret rotation rates 
and light-curve amplitudes together, to obtain information about 
densities and tensile strengths. Instead of analyzing each body 
separately and trying to retrieve information about its density 
and cohesion by using its spin rate and light-curve amplitude 
(and an assumed viewing angle), we address the mean inter- 
nal properties by analyzing the population as a whole. In other 
words, we developed a model that can simultaneously interpret 
the distribution of spin rates and the distribution of light-curve 
amplitudes. 

In this paper, light curve parameters of a vast sample are 
compiled and analyzed. We focus on the analysis of the period 
and shape distribution presented in Sect. 2. Then, with the cur- 
rent data, we complete a correlation analysis with several orbital 
parameters in Sect. 3. Internal properties of the TNOs are derived 
by applying a simple model in Sect. 4, a discussion of the results 
is presented in Sect. 5. Finally, the conclusions are presented in 
Sect. 6. 



2. Period and light-curve amplitude distributions 

One important current limitation is that the light curves of this 
kind of objects are accessible only with medium-size (2 meters) 
or larger telescopes. The telescope time required for this type of 
program is difficult to obtain, and the observed objects are the 
brightest and largest ones. Only several light curves were ob- 
served with the Hubble Space Telescope (HST) reaching a lower 
limit of 20-100 km in diameter ( |Trilling & Bernstein|2006] l. The 
current Kuiper belt object (KBO) absolute magnitude for which 
a good quality light curve can be derived (with a 2 meter tele- 
scope), is H < 9 for TNOs and H< 11 for Centaurs. Smaller 
KBOs can only be observed with space telescopes. 



Using the literature and our recent results (Thirouin et al. 
2009), it is possible to create a database of light curves with ro- 



tational periods and peak-to-peak amplitudes. The special case 
of Pluto, whose spin rate has been altered by tidal dissipation 
caused by Charon, has been removed from the sample. A com- 
pilation of 102 objects, updated on March 2009, is presented in 
Table[T] From that list, only 72 objects have good enough period 
and light-curve amplitude estimates. 

When looking and analyzing the current TNO and Centaur 
light curves, the first noticeable characteristic in most of the de- 
rived light curves is that the amplitudes are very small (<0.15 
mag). The criterium of placing a limit of 0.15 mag was intro- 
duced by Sheppard & Jewitt] (|2002[) and use d by other authors 
(Lacerda & Luu 2006; Ortiz et al. 2003b a) to distinguish the 
"photometrically flat" curves. The brightness changes are prob- 
ably caused by the combined effect of rotation, irregular shape 
and/or albedo markings on the surface. Thus, the abundance of 
low amplitude light curves is indicative of a large abundance 
of quasi-spherical bodies with homogeneous surfaces. However, 
this is to some degree a selection effect because, as mentioned 
above, the true sample is dominated by "bright" objects, which 
are presumably the largest ones and nearly spherical. 

An example of the kind of light curve studied here is pre- 
sented in Fig. [T] In this example, data of the object 2003 VS2 



over several nights was observed and the composite light curve 
was compiled to determine the rotational period and amplitude. 
Unfortunately, it is not always possible to obtain a well deter- 
mined period and/or amplitude. Authors sometimes place a con- 
straint on the values of periods and/or amplitudes, and we note 
that the literature is populated by these kinds of results. In Table 
[T] we list all objects with published rotational period and/or light 
curve amplitude or constraints in one of these values. In any 
case, for all the considerations in the paper only the objects with 
a well determined period/amplitude pair are taken into account. 
That is the case for 72 objects. In case of multiple determinations 
of the period and/or amplitude, the criteria selected the preferred 
value by the author who published that value. If no preferred 
value is mentioned, the mean of the listed values was used. 

2.1. Amplitude distribution 

The KBOs albedo may not be uniform across their surfaces, 
causing the apparent magnitude to vary as the different albedo 
markings on the KBOs surface rotate in and out of our line of 
sight. Albedo or surface variations of the object are usually re- 
sponsible for less than a 20% difference from maximum to min- 



imum brightness of an object (Degew rJ et al.|1979 Magnusson 
|1991| ) in the asteroid case. In Thiro uin et al.| ( |2009| l, it is esti- 
mated that the typical TNO variability caused by albedo mark- 
ings is around 0.1 mag. The reasons for the apparently small 
amount of variability due to albedo markings on the TNO sur- 
faces compared to that of the asteroids may be related to the 
icier surfaces of the TNOs than those of the asteroids. However, 
also the processing of the surfaces is different, because different 
mechanisms can take place. Thus, it is unsurprising that the sur- 
face reflectance variability may be different for TNOs than for 
asteroids. The majority of low-amplitude light curves of large 
objects are believed to be caused by this kind of surface varia- 
tions ( |Sheppard et al.|2008[ ). 

On the other hand, shape variations or elongation of an object 
will cause the effective radius of an object to our line of sight to 
change as the KBO rotates. Smaller KBOs are expected to be 
structurally elongated. To date, few small (20 - 100 km) KBOs 
have been observed with rotational variations, and they appear to 
have larger light curve amplitudes (Trilling & Bernstein 2006). 

In Fig. [2 we show the percentage of objects that have an 
amplitude value within a specific peak-to-peak amplitude (in 0. 1 
mag bin) of the well determined 91 amplitudes (74 for TNOs, 
and 17 for Centaurs) from Table [T] We note that some papers 
only report the amplitude, or some restriction in the amplitude, 
but no period. That is the reason why the number of amplitude 
determinations differ from the rotational period determinations. 

The first noticeable characteristic in that plot is that nearly 
70% of the objects (in both cases, TNOs and Centaurs) have 
light-curve amplitudes smaller than 0.2 mag. The reasons for ob- 
serving a flat light curve could be a spherical body with low con- 
trast albedo marks on the surface, or a pole-on configuration, that 
is, a rotational axis oriented toward the observer. This last option 
is the least probable of all, so it is reasonable to say that most of 
the observed TNOs are "large" round objects with no significant 
inhomogeneities on their surfaces. Objects with flat light curves 
correspond in general to a less collisionally evolved population. 
Of course, the large abundance of the low light-curve ampli- 
tude population could be to some degree an observational bias 
because the KBOs with known rotational properties are mainly 
the brightest ones. For the Centaurs, the vast majority of objects 
have light-curve amplitudes between 0.1 mag and 0.2mag, which 
is indicative of some elongation in the shape of the bodies. 
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2.2. Period distribution 

In the following analysis, all objects with known rotational pe- 
riod are considered. An exception was made for Pluto, which 
was excluded. This exception is because the long rotational pe- 
riod (due to the tidal locking of the system) cannot be represen- 
tative of the rotational properties of the sample. 

From KBO light curves alone it is difficult or even impos- 
sible to determine whether the variability is caused by either 
albedo variations or an elongated shape. Therefore, the rotation 
periods derived from the photometric periods often have an un- 
certainty of a factor 2. A double-peaked periodic light curve is 
expected to be seen in a non-spherical body, since the projected 
cross-section should have two minima (short axis) and two max- 
ima (long axis) during one complete rotation of the KBO. Light 
curves caused by surface inhomogeneities are single peaked. 
Figure [3] (in black) shows the histogram of rotational frequen- 
cies (single-peak rotational period) for all the objects found in 
the literature. This histogram was compiled by assuming that all 
objects have single-peaked light curves, which is obviously not 
true and is an extreme case. For the well known cases in which 
the light curves are certainly double-peaked, we used the rotation 
period corresponding to half the true rotation period. Therefore, 
the resulting distribution of rotation periods has a short mean 
value, which is a crude lower limit to the mean rotation period. 
A Maxwellian curve was fitted and overplotted in the histogram 
for a mean period of 5.13 hours. 

The obvious next exercise is to plot the histogram of rota- 
tional frequencies for all objects with double-peak light curves 
(in other words, we assume that all rotation periods are twice the 
periods derived by a basic periodogram analysis). This is again 
an unrealistic result, but provides an upper limit to the rotational 
period distribution. This is the same as considering all the light 
curves to be produced by non-spherical bodies. The Maxwellian 
fit infers a mean frequency of 2.47 cycles/day or 9.70 hours (data 
in grey). 

Finally, a more reliable histogram (presented in Fig. |4]i was 
compiled by assuming that all the light curves with amplitudes 
smaller than or equal to 0. 15mag are single peaked (this is equiv- 
alent to assuming that the main cause for their rotational varia- 
tion is albedo markings) and those with amplitudes larger than 
0. 1 5mag are double peaked. Thus, with this criterium we remove 
the factor of 2 ambiguity and derive a more realistic distribution. 
There is also the possibility of deciding on a case-by-case basis, 
but some cases are ambiguous. 

If we calculate the raw average rotational period, the ob- 
tained values from Table[T]are 6.95, 6.88, and 6.75 hours for the 
entire sample, for both the TNOs and the Centaurs only, respec- 
tively. On the other hand, if we fit Maxwellian distributions to 
the rotational frequencies, the mean of the distributions change 
slightly. The two main populations are separated here, giving a 
mean rotational period of 7.35 hours (3.19 cycles/day) for the 
entire sample (upper panel), 7.71 hours (3.11 cycles/day) for the 
TNOs (center panel), and 8.95 hours (2.68 cycles/day) for the 
smaller Centaurs (lower panel). 

As we can see in Fig. |4] there appears to be some fast ro- 
tators. In the center panel, for the TNOs, there is Ceto, 1996 
TP66, and 2004 GVy with single peak rotational periods of 2.21, 
1.96, and 2.93 hours, which are quite short and unrealistic. The 
reason for this is that all 3 of the afore mentioned objects have 
light-curve amplitudes smaller than 0.15mag and were consid- 
ered to have single-peaked light curves, whereas in reality they 
have double-peaked light curves. A similar case occurs for the 
Centaur Chiron and 2002 GZ32 (lower panel), which have a sin- 



gle peak rotational period of 2.96 and 2.9 hours, respectively, for 
the same reason. 



3. Correlations 

To search for correlations between different physical and orbital 
parameters (rotational period, peak-to-peak amplitude, semi- 
major axis, perihelion distance, aphelion distance, eccentricity, 
inclination), we use non-parametric statistical tests. 

The presence of several periods and/or amplitudes for most 
of the objects, encouraged us to develop a customized code that 
randomly chooses one combination. The program created the in- 
dividual combinations of period and amplitude pairs and calcu- 
lated the correlation between one of these combinations (of pe- 
riod and amplitude) and the selected variable (e.g., colors, orbital 
parameters). The final choice of correlation, presented in Table 
[2] is the median among those correlations. 

The statistical methods do not assume any particular popu- 
lation probability distribution, nor any functional shape. In par- 
ticular, we compute all the correlations using the Spearman rank 
correlation, p ( Spearman 1 1904] >. This method is distribution-free 
and less sensitive to outliers than others methods. We studied 
the strength of the correlations by computing the Spearman co- 
efficient p and the significance level (SL). The p coefficient has 
values between -1 and 1. A correlation may exists if p > 0, while 
p < implies a possible anti-correlation. A zero value means 
that there is no correlation at all. The significance of p is mea- 
sured by SL, the probability that the null hypothesis (samples 
not correlated) is not true. We consider that correlations or anti- 
correlations are strong when p (in absolute value) is greater than 
0.6, a possible weak correlation exists when p values are be- 
tween 0.3 and 0.6, and no correlations are present for values less 
than 0.3. The following criteria were used: SL greater than 95% 
= reasonably strong evidence of correlation; SL greater than 
97.5% = strong evidence of correlation; and SL greater than 99% 
= very strong evidence of correlation. The interpretation of p, 
SL, and the graphical representation of the magnitudes allows us 
to detect possible correlations among the different magnitudes. 
Apart from these criteria, we consider that a possible correla- 
tion could exist when p is greater than or equal to 0.3 and SL 
is greater than 80%, in particular for samples of only a few ob- 
jects. We emphasize that this particular threshold (p >=0.3; SL> 
80%) provides only a hint of a possible correlation. We note that 
in the majority of these cases, the sample is so small, that more 
observations are required in some of the cases listed in Table [2] 

The Spearman statistical analysis of the results was made by 
separating the dataset into different dynamical subgroups, i.e., 
classical objects, which exhibit orbits with moderate eccentrici- 
ties; resonant objects, trapped by Neptune in mean motion res- 
onances (the Plutinos being the most representative population, 
in the 3:2 resonance), and scattered disk objects (with orbits of 
high eccentricities and sometimes high inclinations, due to close 
encounters with the planet Neptune in the past). The classical 
group is sub-classified into two sub-groups: the so-called (dy- 
namically) hot classical objects with orbital inclinations > 4.5°, 
and the (dynamically) cold classical objects, with inclinations < 
4.5°. Finally, the Centaurs represent TNOs scattered towards the 
inner Solar System and they reside between the orbits of Jupiter 
and Neptune. 

The analyzed parameters were the rotational period, light- 
curve amplitude, colors (B-V, R-I and V-R), absolute magnitude 
H, and orbital parameters (a,e,i,Q and q). It should be men- 
tioned that a data set of colors was presented in |Santos-Sanz| 
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|et al.| l |2009) . All parameters were tested against all the other pa- 
rameters in the entire sample and then in sub-populations. The 
corresponding correlation coefficients and significance levels are 
listed in Table [2] where only the relevant cases are shown. We 
made the test correlations with all the parameters in all the pop- 
ulations and then listed only the correlations that were found to 
be significant. 

For the case of amplitude versus H, values of p and SL for 
all the sample and in particular for the classical TNOs, clearly 
indicate that the smaller (and collisionally evolved) objects are 
probably more elongated than the larger ones. This could explain 
the larger amplitudes observed for smaller objects. In the classi- 
cal population (39 objects), the value of p = 0.47 and SL = 99.5 
% indicates a stronger correlation than the entire population with 
p - 0.40 and SL = 99.8 %. For the Centaurs and Plutinos, we do 
not find any evident correlation between these two magnitudes. 

The particular case of P versus B-V correlations may sug- 
gest that objects with shorter rotation periods may have suffered 
more collisions (and thus have bluer surfaces) than objects with 
longer ones. The collisions could resurface the objects with fresh 
ices making them bluer and modifying the spins toward faster ro- 
tation periods. There is a strong evidence of a correlation in the 
classical subpopulation, although this sub-sample has fewer than 
20 objects. More color determinations of objects with known pe- 
riods would clearly be helpful, not only to confirm or exclude 
this possible correlation but to confirm whether there is some 
color variation in the surface of the body. The results obtained for 
the P versus H relations of Centaurs (correlation) and Plutinos 
(anti-correlation) might be related to the different mean sizes of 
the Centaur and Plutino populations and could be related to a 
different primordial origin and/or evolution. In particular, phys- 
ical parameters such as the color B-V and rotational period are 
known only for brighter and larger objects. We, again note the 
low number of objects. We found some weak anti-correlations 
between the amplitudes and orbital parameters Q and e, whose 
possible causes are not clear to us. 

Finally, another interesting plot is the rotational period ver- 
sus absolute magnitude (an indicator of the size), where we can 
see (Fig. [5]l that a rotation barrier could exist close to 4 hours. 
All the points have an error bar that indicates the period deter- 
mination error when the result is unique and a larger error when 
the period has several values in the literature. A typical example 
of multiple determinations in the literature is Eris, for which the 
period has been determined to be either 13.7 h (Duffard et al. 
[20081 or 25.92 h ( |Roe et al.|2008] >. 

4. Densities and internal structures 

With the rotational period sample, it is possible to derive some 
constraints on the density of the objects. Several authors have 
presented different estimations on TNO densities. Pravec & 
Harris (2000) derived a simple expression that approximates the 
critical (minimum) period for a rotating body as a function of 
density and light-curve amplitude. This relation assumes a fluid 
body. Sheppar d et al.| ([2008) used the periods and shapes (de- 
rived from light curves) to identify an apparent trend of larger 
(brighter) objects being denser. Finally, Tancredi & Favre (2008 ) 
analyzed TNO densities and internal structures to determine 
whether TNOs could be dwarf planets. 

We note that objects with radii r > 200 km have presum- 
ably never been disrupted by impacts (Davis & Fari nella|1997| l, 
but have probably been in some way fractured. If we assume 
that TNOs have relaxed to equilibrium shapes, then their rota- 
tion states can be used to place limits on their densities. The 



centripetal acceleration produced by self-gravity (bulk density) 
must be sufficient to hold the material together against the in- 
ertial acceleration due to rotation. Although TNOs consist of 
solid material, their presumably fragmentary structure validates 
the fluid approximation as a limiting case. 

Densities have been proposed to be low, 600 kg/m 3 for 
2001 QG298 and in some cases, covering a wider range to higher 
densities of 2600 kg/m 3 as in 2003 EL 61 . The lower density is 
similar to that inferred for some comets, which are objects two 
orders of magnitude smaller. Bulk densities substantially below 
that of pure water ice are usually taken to imply an ice-rich com- 
position and porosity. It is expected that density variations are 
not only related to intrinsic compositional differences, but might 
be caused by variations in the porosity. On the other hand, for the 
dwarf planet class, porosity cannot be important. It needs to be 
mentioned that a density of 2600 kg/m 3 implies a rock/ice ratio 
of ~85/15% when the internal structure is modelled (McK rnnon| 
|et al.|2008| >. 

4.1. A simple Monte Carlo model 

Taking into account all the light curve amplitude and rotational 
period data presented, it can be seen that "flat curves" are dom- 
inant. As mentioned before, this could be because we are ob- 
serving large presumably round objects. In the near future, more 
light curves should be obtained of smaller objects and the ques- 
tion may then arise of whether this tendency remains and how 
many non-flat curves are observed. To estimate how many non- 
flat curves are expected in an observed population, we developed 
a simple model. 

The first step in our model is to assume a Maxwellian rota- 
tional frequency distribution such as that obtained in this work 
(Fig. |4|, from the observational data. Then, each object from a 
set of 100 000 of fixed density is randomly assigned to a rota- 
tional period from the distribution. The Maxwellian rotational 
frequency distribution adopted ranges from 8 to 0.5 cycles/day, 
with a bin of 0.1 and a peak at 3.11 cycles/day (7.71 hours), in 
the observed distribution. 

One important constraint of the model is that light curves 
are assumed to be related to the shape of the body. The model 
aims to determine the percentage of objects that create non-flat 
light curves due to elongated shapes. Other causes of brightness 
variation (not assumed here) are, for example, eclipses in binary 
objects or surface albedo variations. 

For the model calculations, it is assumed that each object 
rotates in hydrostatic equilibrium and the Chandrasekhar for- 
malism (Chandrasekhar 1987) is used. The questions that we try 
to answer is how many objects (of a given density) with differ- 
ent rotational periods could have an irregular shape with axes of 
lengths such that a>b>c (a Jacobi ellipsoid), and how many ob- 
jects acquire a MacLaurin shape (i.e., a=b>c). The final percent- 
age of MacLaurin or Jacobi objects should be compared with the 
observed data. 

Jacobi ellipsoids produce a non-flat light curve if the ratio of 
the axis lengths are greater than the photometric error. It is then 
necessary that the rotating axes are not aligned along the ob- 
server line of sight. MacLaurin objects will produce a flat light 
curve independently of the rotational axis direction. In all cases 
the objects are considered at opposition and in a simple rota- 



tion around the shortest axis. According to Burns & Safronov 
( |1973| l, rotational excitation damping times are orders of magni- 
tude shorter than the age of the Solar System and are also orders 
of magnitude shorter than the typical time between collisions in 
the transneptunian region. In their study of ellipsoids, |Molina| 
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|et al.| ( [2003] l emphasized that damping times can be even shorter 
if one takes into account the energy dissipation caused by inter- 
nal stresses. Thus, it is expected that the transneptunian objects 
must be purely in rotation around their principal axes of max- 
imal moment of inertia, because the collisions that can excite 
their rotations take far longer to occur than the typical rotational 
excitation damping time. The TNOs have been found to be in 
pure rotation in other works, e.g., Ortiz et al. (2003a); Lacerda 
|& Luu| ( |2003l ) 

For all objects in each run, the dimensionless parameters O 
and A (dimensionless angular velocity and dimensionless an- 
gular momentum) ( Chandrasekha r|T987 1 were calculated. From 
these pairs of parameters and using Chandrasekhar's tables, it 
is possible to determine whether the object has a Jacobi or 
MacLaurin figure or if a non-equilibrium figure is possible. To 
obtain a complete sample of all possible densities, the model was 
run for several density values between 500 and 3000 kg/m 3 . 

4.2. Results for the model 

For each density, the percentage of Jacobi, MacLaurin, and non- 
equilibrium figures were calculated and are presented in Table 
[3] In the case of the Jacobi ellipsoids, there is an interval of al- 
lowed rotational periods where the objects are in equilibrium. 
For example, for a density of 1000 kg/m 3 , only 12.61% of the 
objects are Jacobi ellipsoids and they rotate with periods be- 
tween 6.31 and 7.05 hours. For this same density, 55.63% of 
objects have a MacLaurin figure and the remaining 31.75% are 
non-equilibrium figures that were discarded from the sample. 

In Fig. [6] the percentage of Jacobi, MacLaurin, and non- 
equilibrium objects is plotted against density for all the runs. 
As can be seen, when the density increases the percentage of 
MacLaurin spheroids and non-equilibrium figures, increase and 
decrease, respectively. The percentage of Jacobi ellipsoids has a 
maximum close to the densities of 1200-1300 kg/m 3 . 

The amplitude of the light curve of each synthetic object is 
then calculated from the axis ratio b/a of the Jacobi ellipsoid. 
For a fixed density, the ratio b/a of the ellipsoid increases with 
the rotational period. The obtained axis ratio was then modified 
by a random orientation of the rotational axis. This means that an 
elongated Jacobi ellipsoid would be seen pole-on and no light- 
curve variation would be observable. In Fig. [7] the light curve 
amplitude for several runs of the model is plotted against rota- 
tional period. As can be seen, for each density there is an in- 
terval of allowed rotational periods for Jacobi ellipsoids. All the 
period-lightcurve amplitude pairs below that curve are allowed 
after the correction of the viewing angle. We note that the al- 
lowed rotational period interval for a lower density is longer 
than that allowed for higher densities, as expected. The rotational 
period interval becomes shorter when the density increases, in 
other words, the curve becomes steeper as the density increases. 

Observed data is included in the figure to compare with the 
modelled data. From this figure, the density of a real object can 
be estimated by assuming that the objects are viewed equatori- 
ally. Haumea has a estimated density of 2500 kg/m 3 , while most 
of the other selected objects have variations between densities 
500 to 1000 kg/m 3 . 

Finally, in Fig.[8]the percentage of objects with different fig- 
ures is presented for several densities. MacLaurin figures were 
added to the Jacobi ellipsoids in the case of small light curve 
amplitudes (0.1 mag). Here it is assumed that if a light curve 
amplitude is detectable in a MacLaurin object, it will be smaller 
than ~ 0.1 mag, because it will be produced by albedo features. 
For the different model runs (densities) we can observe that the 



percentage of objects with small light curve amplitude (O.lmag) 
is always higher than the percentage of the more elongated ones, 
as expected. This percentage increases from 30% for a low den- 
sity such as 500 kg/m 3 to more than 90% for a higher density like 
2500 kg/m 3 . On the other hand, as can be seen, the percentage 
of objects with large amplitudes is independent of density and is 
always smaller than 20%. The percentage of MacLaurin objects 
sharply increases with density and this is the main reason for the 
differences between the curves at the left part of Fig. [8] (small 
light curve amplitudes). 

In Fig. [8] the observational data are also plotted separated 
into two intervals of absolute magnitudes. The percentage of ob- 
jects in the interval of H between [-1,5] and [-1,7] is taken from 
Table Q] The behavior of both curves are the same as for the 
modelled ones. A large percentage of the objects have low am- 
plitude light curves. Comparing this behavior with the modelled 
curves it is possible to estimate the mean density of the entire 
sample. As can be seen, the observational data agree with the 
output of the corresponding model run for 1000 kg/m 3 density. 
When smaller objects are included (black squares), the observed 
points are reproduced most accurately with an even lower den- 
sity. 

At the other side of the plot, it can be seen that observational 
data are lacking with respect to the model prediction. A blow- 
up of the plot for the larger amplitudes is inserted in Fig. [8] As 
can be seen, some objects with larger light curve amplitudes are 
expected to be found. 

To summarize, if hydrostatic equilibrium is considered with 
a Maxwellian rotational frequency distribution, a large percent- 
age of "flat light curves" is obtained. This is because of the large 
percentage of MacLaurin figures added to the Jacobi ellipsoid 
with small amplitude (due to a small axis ratio or viewing an- 
gle). It is important to note that for the results mentioned above, 
only a distribution of rotational frequencies and a fixed density 
are assumed. No other physical variable is needed. 

Finally, if a size distribution is included in the model, it is 
possible to compare the objects with the real sizes. For that, a 
cumulative surface density of TNOs was included as described 
in |Petit et al.| ([2008 ) between 18 < R < 22 magnitude. Assuming 
a mean albedo of 0. 12 and mean distance of 40 AU, the smallest 
objects in this group have a radius of 86 km. In other words, hy- 
drostatic equilibrium seems to be valid for objects smaller than 
assumed before ( Shepp ard et al.| [2008 ). This could be because 
the internal cohesion forces are smaller in TNOs than in aster- 
oids. An estimated value of the tensile strength for this size that 
is overcome by the object's gravity, allowing the objects to be- 



have like a fluid, is obtained from equation 3 of Thirouin et al. 
(2009). The tensile strength is only 0.1 MPa. This small tensile 
strength is much smaller than that of typical geophysical solids. 



5. Discussion 

With all the available data, we have compiled different his- 
tograms to see the distribution of physical parameters such as 
rotational periods and light-curve amplitudes. With this larger 
number of light curves, the results obtained for the mean ro- 
tational period and peak-to-peak amplitudes are more reliable 
than in previous studies. The criteria to separate flat light curves 
from the others is important when selecting a single or double 
peak rotational period. As is demonstrated, the distribution of 
periods infers a lower limit (for a single-peak period) and upper 
limit (for double-peak periods). The result of the combination of 
both periods is a more reliable result. The threshold of 0.15mag 
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for the transition from the albedo-marking light curves to shape- 
induced light curves is somewhat arbitrary and all the results 
are based on this value. If a light curve amplitude magnitude of 
O.lOmag is considered the obtained results are similar, reaching 
a mean of 7.84 hours and 7.12 hours for all the sample and for 
only TNOs, respectively. The goodness of the Maxwellian fit is 
similar to the 0.15mag case. However, if we use a threshold of 
0.2mag, the fit is degraded. 

Why are their surfaces so uniform? The sample of flat light 
curves is dominated by larger objects. These bodies are expected 
to be round and in hydrostatic equilibrium. In these bodies, it 
is natural to ask: could the reason for flat light curves be the 
presence of a thin atmosphere or recent outgassing/resurfacing? 

The environment in the trans-neptunian region might favor 
the formation of tenuous permanent or transient atmospheres. 
This idea is supported by several pieces of evidence such as 
(i) the intrinsic and/or induced activity on some bodies, such 
as in Chiron ( Duffard et al.| |2002) and other Centaurs. Activity 
has been reported for comets at large heliocentric distances and 
much beyond the water ice sublimation limit ( |Hainaut et al. 
2000 ). (ii) Abundance of volatiles with sublimation temperatures 
comparable to the local temperature, (iii) The most important, 
the atmosp heres of Pluto ( |Yelle & Eliiot|1997[ > and Triton ( jYelle 
|et al.|1995| l. 

Our study is based on a photometric analysis, so we cannot 
determine the presence of atmospheres. The observational result 
is that most TNOs have small amplitudes. Taking the evidence 
together, it appears reasonable to expect the formation of thin 
gaseous and/or dusty atmospheres in large TNOs, because they 
would have the mass needed to avoid direct escape of materi- 
als into space. There is a study by Lykawka & Mukai ( 2005 ) 
about the theoretical presence of an atmosphere according to 
several criteria and the possibility of the presence of a thin at- 
mosphere in the studied bodies that takes into account volatiles 
such as methane, nitrogen, and carbon monoxide, the most com- 
mon chemical species in the outer Solar System. It is noticeable 
that only TNOs with a diameter of 700 km (albedo criteria are 
explained in that article) can retain an atmosphere. The large 
self-gravity of the TNOs should allow them to remain close to 
hydrostatic equilibrium, retain extremely volatile ices, possibly 
have tenuous atmospheres, and be differentiated. Thus, the sur- 
faces as well as the interior physical characteristics of the largest 
TNOs may differ significantly from those of the smaller objects. 

A nearly uniform surface may be explained by an atmo- 
sphere that is frozen onto the surface of the objects when close 
to aphelion. The atmosphere, such as Pluto's, may become ac- 
tive when near perihelion, effectively resurfacing itself every few 
hundred years. Taking into account all the known light curves 
with known amplitudes and the work of Lykawka & Mukai 
(2005), it is possible to identify possible candidates to be spec- 
troscopically observed by searching for spectral signatures of an 
atmosphere. 

The transient atmospheres and atmospheric freeze-out as a 
means of resurfacing and homogenization would apply to the 
largest objects only. An alternative scenario might be the colli- 
sional resurfacing, which may explain other phenomena in the 
transneptunian belt ( |Gil-Hutton|2002l [GiFHutton et al.|2009] >. 

In the model presented here, for different runs of only one 
density, the percentage of MacLaurin figures is always superior 
to 50% for densities higher than 900 kg/m 3 . As can be seen in 
Fig. [8] the percentage of observed objects with light curve am- 
plitudes smaller than 0.2mag is always 70%. When compared to 
real data, it needs to be mentioned here that this sample is biased 



as observers tend to publish light curves with large amplitudes 
and easy-to-measure rotational periods. 

According to our simple model, lower mean densities are 
obtained when we analyze the sample with the largest size span, 
from H=[-l,7], compared to the interval H=[-l,5]. Hence, it ap- 
pears that the "small" objects are also in hydrostatic equilibrium, 
but have lower densities than the largest objects. Therefore, an 
obvious refinement of our model could be to incorporate a size- 
dependent density together with the size distribution of the bod- 
ies whose brightness variability has been studied. Other im- 
provements in the model could be the use of non-random spin 
axis orientations to test the effect of different spin axis distribu- 
tions. 

Despite the simplicity of our model, its results are compat- 
ible with observations. Currently, the only input of the model 
is a Maxwellian rotational frequency distribution, but the real 
distribution might differ from a Maxwellian, since this func- 
tion is representative of a collisional equilibrium population. 
However, the true scenario might be represented by a combina- 
tion of a collisional evolved populations and a more primordial 
one. Therefore, an additional improvement of the model would 
be the possibility to include different rotational frequency distri- 
butions. 

As a caveat to the densities derived here, we can mention the 
work of Holsapple (2001 ) and Holsapple (2004) who determined 
the spin limits of bodies using a model for solid bodies without 
tensile or cohesive strength. That theory included the classical 
analysis of fluid bodies given by MacLaurin, Jacobi, and others 
as a special case. For the general solid bodies, it was shown that 
there exists a very wide range of permissible shapes and spin 
limits, and explicit algebraic results for those limits were given. 
Holsapple (2007) extended these analyse to include geological- 
like materials that also have tensile and cohesive strength. In that 
paper, the author asked whether larger bodies such as the TNOs 
should have relaxed fluid shapes or might have shapes that are 
allowable by granular zero cohesion material? The important re- 
sult is that if they are not governed by fluid theory, the mass and 
shapes that have been inferred could be substantially incorrect 
Polsapple|2007] ). 

Finally, we must mention that the search for correlations be- 
tween light curve variables and orbital parameters is hampered 



by bias and selection effects ( Binzel et al. 1989 1, and our conclu- 
sions in this respect are weak. However, the correlations between 
light curve parameters and physical properties such as colors and 
size that we derived are less affected by bias and we believe are 
reliable. 

6. Conclusions 



We have analyzed the light-curve data of Thirouin et al. ( 2009 ) 
and others in the literature to derive light-curve parameters, 
which have been compiled here. We have presented the most 
complete distribution of rotation periods and amplitudes so 
far. For the rotation periods, the mean value obtained by a 
Maxwellian fit in rotational frequencies is 7.71 hours (3.11 
cycles/day) considering only the TNOs, 8.95 hours (2.68 cy- 
cles/day) considering only Centaurs, and 7.35 hours (3.19 cy- 
cles/day) considering all the sample altogether. These results 
were obtained by taking into account the criteria of consider- 
ing a single peak light curve for those objects with amplitudes 
lower than 0.15 mag and a double peak light curve for the more 
elongated ones. 

This mean rotational period is slightly longer than the mean 
value for the largest Main Belt asteroids, 6.9 hours ( |Binzel et al.| 
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1989), but we state that this value may be affected by observa- 
tional biases. It is important to stress that in photometry devoted 
to the study of the rotational properties of minor bodies, there 
exists and important observational bias towards shorter periods 
and higher amplitude light curves. How this bias can be esti- 
mated remains an open problem, but it certainly depends on the 
number of objects for which a period is obtained. Some steps 
in this direction were already taken by Thirouin et al. (2009). 
Concerning other biases, if observational time for larger tele- 
scopes is granted then smaller TNOs will be accessible for the 
determination of the light curves, and the sample will not only 
increase in size but the bias in size will become less important. 

One of the most important correlations found in this work 
is the light curve amplitude and absolute magnitude of all the 
sample of KBOs. This correlation is not a surprise and is prob- 
ably related to the most collisionally evolved small population, 
which should be more elongated objects. On the other hand, the 
less collisionally evolved population (larger and then brighter 
objects) is the more numerous one in our sample (~60%), so we 
caution the reader that some of our results could be somewhat bi- 
ased. Another interesting correlation is the rotational period and 
B-V color, which could indicate the age of the surface (younger 
surfaces being bluer). This correlation appears to be consistent 
with the idea that the more collisionally evolved objects spin 
faster and have bluer surfaces. 

A simple shape model has been developed that is able to re- 
produce some results on the relative percentage of presence of 
MacLaurin/Jacobian figures, when assuming hydrostatic equi- 
librium. In the entire sample, less than 12% of the objects in 
hydrostatic equilibrium have a Jacobi shape. Most of the objects 
have MacLaurin figures that have small light curve amplitudes 
(because of albedo marks). The model that best fits the observa- 
tions requires densities between 1000 to 1200 kg/m 3 . We must 
mention that as was highlighted in Thiroui n" et al.| ([2009) there 
is an observational bias in the percentage of low amplitude light 
curves, which means that perhaps 75% of the objects have low 
light curve amplitudes. This implies that the mean density is 
likely somewhat higher, around 1500 kg/m 3 (see Fig. [8]). 

For plausible albedos of 0.04 to 0.20, the absolute magnitude 
range H=[-l,7] includes objects with diameters as small as ~ 120 
km. Those objects would also be in hydrostatic equilibrium. We 
end with the provocative thought that the great majority of ob- 
served TNOs would qualify to be dwarf planets, because they 
appear to meet the hydrostatic equilibrium condition. 
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Fig. 1. Example of a KBO light curve from 
obtained during several nights. 
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Fig. 2. Percentage of objects vs. peak-to-peak light curve ampli- 
tude for all TNOs (black) and Centaurs (grey). 
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Fig. 3. Histogram in rotational frequencies for all the objects 
considering only a simple period rotation (black) and double pe- 
riod rotation (grey). A maxwellian fit to the single peak period 
for all the object gives a mean not-reliable lower period of 5.13 
hours (black bars/line). Double peak period is also represented 
(grey bars). The maxwellian fit to the data 
gives a mean not-reliable upper period of 9.70 hours. 
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Only Centaurs, simple/double, n=1 4 
Mean = 2.68 cycles/day = 8.95 hours 
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Fig. 4. Histogram in cycles/day for the whole sample (top panel), 
the TNOs (center panel) and Centaurs (lower panel). From 
maxwellian fits to the rotational frequencies distribution the 
mean rotation rates are 7.35 hrs for the entire sample, 7.71 hrs 
for the TNOs alone and 8.95 hrs for the Centaurs. We chose to 
use the single peak period for the objects whose amplitude is 
<=0.15 mag and double peak period for the objects whose am- 
plitude is >0.15 mag. 
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Fig. 5. Rotational Period against absolute magnitude for all the 
objects. The horizontal line at 4 hrs is the spin barrier mentioned 
in the text. 
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3000 3500 



4000 



Duffard, R. et al.: KBOs CCD photometry 




Fig. 7. light curve amplitude vs. rotational period for several 
model runs compared with some observational data. Densities 
are in kg/m 3 
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Fig. 8. Percentage of the objects having a certain light curve am- 
plitude from the model for different densities. MacLaurin object 
and Jacobi ellipsoid were added together only for a light curve 
amplitude of 0.1 mag. Circles and squares are observational data 
using values of absolute magnitude between [-1,5] and [-1,7], 
respectively. 
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Table 1. List of TNOs whose variability has been studied. When a multiple determination of the rotational period is present, the preferred period 
is marked in bold. 



Object 


Designation 


Single Peak 


Double Peak 


Amplitude 


Absolute 


Ref b 






Period [h] 


Period [h] 


[mag.] 


magnitude 11 




(134340) Pluto 




153.6 




0.33 


-0.7 


B97.TH90 


Charon 




153.6 




0.08 


0.9 


B97 


(148780) Altjira 


2001 UOiR 






<0.10 


5.6 


S07 


(55576) Amycus 


2002 GB 10 


9.76 


19.52 


0.16±0.01 


7.8 


Th09 


(8405) Asbolus 


1995 GO 


4.47 


8.9351±0.003 


0.55 


9.0 


D98a,K00 










0.34 




RT99 


(54598) Bienor 


2000 QC243 


4.57±0.02 


9.14±0.04 


0.75±0.09 


7.6 


O03b 






9.17 




0.34±0.08 




R07a 


(66652) Borasisi 


1999 RZ 253 






<0.05 


5.9 


LL06 


(65489) Ceto 


2003 FX 12 s 




4.43±0.03 


0.13±0.02 


6.3 


SJ02 


(19521) Chaos 


1998 WH24 






<0.10 


4.9 


SJ02.LL06 


(10199) Chariklo 


1997 CU26 






<0.1 


6.4 


LL06 


(2060) Chiron 


1977 UB 




5.9180 ± 0.0001 


0.088±0.003 


6.5 


B89 


(83982) Crantor 


2002 GO, 


(6.97 or 9.67)±0.04 




0.14±0.04 


9.1 


O03b 










0.34 




RT99 


(60558) Echeclus 


2000 ECc, g 


13.401 


26.802 


0.24±0.06 


9.0 


R03 


(31824) Elatus 


1999 UG 5 


13.25 




0.24 


10.1 


G01 






13.41±0.04 




0.102±0.005 




B02 


(136199) Eris 


2003 UB 3n 


13.69/28.08/32.13 




<0.1±0.01 


-1.2 


D07 






3.55? 




-0.05 




L07 










<0.01 




R07a,S07 






25.92 




0.1 




R08 


(136108) Haumea 


2003 EL(,i 




3.9154±0.0002 


0.28 ±0.04 


0.2 


R06 








3.9155±0.0001 


0.29±0.02 




LJ08 






1.96 


3.92 


0.23 ±0.02 




Th09 


(10370) Hylonome 


1995 DW 2 






<0.04 


8.0 


RT99 


(38628) Huya 


2000 EB173 


6.68/6.75/6.82 




<0.1 


4.7 


O03b 










<0.06 




O03b,SJ02,S02 










<0.04 




SJ03.LL06 


(28978) Ixion 


2001 KX 76 






<0.05 


3.2 


O03b,SJ03 










<0.10 




SJ02.LL06 


(136472) Makemake 


2005 FYg 


11.24±0.01 


20.54/22.48 


0.03±0.01 


-0.3 


O07 






7.65 


15.30 


0.014±0.002 




Th09 


(7066) Nessus 


1993 HA 2 






0.5 


9.6 


RT99 










<0.2 




D98b 


(52872) Okyrhoe 


1998 SG 15 


4.86/6.08 


9.72/12.16 


0.07 ±0.01 


11.3 


Th09 


(90482) Orcus 


2004 DW 


7.09/10.47±0.01/17.43 


20.94 


0.03±0.01 


2.3 


O06 






13.19 




0.18 




R07a 










<0.03 




S07 






10.47 


20.94 


0.03±0.01 




Th09 


(5145) Pholus 


1992 AD 




9.98 


0.15/0.60 


7.0 


B92,H92,F01,T05 










0.15 




RT99 


(50000) Quaoar 


2002 LM 60 




17.6788±0.0004 


0.13±0.03 


2.6 


O03b 






8.84 




0.18±0.10 




R07a 






9.42 


18.84 


-0.3 




L07 






8.84 


17.68 


0.13±0.04 




Th09 


(90377) Sedna 


2003 VB 12 


10.273±0.003 




0.02 


1.6 


G05 


(8861 1) Teharonhiawako 


2001 QT 297 






<0.15 


5.5 


Os03 


(32532) Thereus 


2001 PTn 


4.1546±0.0001 


8.3091±0.0001 


0.16±0.02 


9.0 


O03b 








8.34 


0.16 




FD03 








8.34 


0.34±0.08 




R07a 








8.4 


0.15 




F01 


(42355) Typhon 


2002 CR 46 


9.67 


19.34 


0.06±0.01 


7.2 


O03b 










<0.05 




SJ03 






>5 








D08 






9.67 


19.34 


0.06±0.01 




Th09 


(20000) Varuna 


2000 WR 106 


3.1718 


6.3436±0.0001 


0.43±0.01 


3.6 


O03b 








6.34 


0.42±0.03 




SJ02 






3.17 


6.34 


0.5 




F01 








6.34358±0.00002 


0.42 




B06 








6.344 


0.49±0.17 




R07a 






3.1709 


6.3418 


0.43±0.01 




Th09 


(15789) 


1993 SC 


7.7 




0.04 


6.9 


T97 










<0.04 




RT99 


(15820) 


1994 TB 


3.0/3.5 


6.0/7.0 


0.26/0.34 


7.1 


RT99 










<0.04 




SJ02 


(19255) 


1994 VK 8 


3.9/4.3/4.7/5.2 


7.8/8.6/9.4/10.4 


0.42 


7.0 


RT99 






4.75 








CB99 


(24835) 


1995 SM55 


4.04±0.03 


8.08±0.03 


0.19±0.05 


4.8 


SJ02 


(32929) 


1995 QY 9 


3.5 


7.0 


0.60 


7.5 


RT99 






Between 3.3 and 3.7 




0.60 




RT99 








7.3±0.1 


0.60±0.04 




RT99.SJ02 


(26181) 


1996 GQ 21 






<0.10 


5.2 


SJ02 


(15874) 


1996 TL 66 


12.1 




<0.12 


5.4 


O06 










<0.06 




LJ98.RT99 






6/8.04/ 12 


24 


0.09±0.02 




Th09 


(19308) 


1996 T0 6( , 


3.96±0.04 


7.92±0.04/5.9/9.6 


0.25±0.05 


4.5 


SJ03 








11.9 


0.25±0.05 




B06 








6.25±0.03 


0.12 to 0.33 




H00 










<0.10 




RT99 


(15875) 


1996 TP«, 


1.96 




<0.04 


6.8 


CB99 










<0.12 




RT99 
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Table 1. continued. 
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Object 


Designation 


Single Peak 
Period [h] 


Double Peak 
Period [h] 


Amplitude 
[mag.] 


Absolute 
magnitude' 1 


Ref b 


(118228) 


1996 TQ 66 


- 


- 


<0.22 


7.1 


RT99 




1996 TS 66 


- 


- 


<0.14 


6.5 


LL06 






_ 


_ 


<0.16 




RT99 


(79360) 


1997 CS 2 , 




- 


<0.08 


5.1 


SJ02 








- 


<0.22 




RT99 




1997 CV 29 


- 


15.8 


0.4 


7.3 


CK04 


(33128) 


1998 BU 4g 


(4.9 or6.3)±0.1 


(9.8 or 12.6)±0.1 


0.68±0.04 


7.2 


SJ02 


(91133) 


1998 HK 15 i 


- 


- 


<0.15 


7.6 


SJ02 


(26308) 


1998 SMi 65 


- 


7.1±0.1 


0.45±0.03 


5.8 


SJ02 






3.983 


- 


0.56 




R01 


(35671) 


1998 SN 165 


- 


8.84 


0.16±0.01 


5.8 


LL06 






5.03 


- 


0.15 




P02 


(33340) 


1998 VG 44 


- 




<0.10 


6.5 


SJ02 




1998 XY, 5 


- 


- 


<0.1 


6.2 


CB01 






1.31 


- 


-0.1 




CB01 


(26375) 


1999 DE 9 


>12? 


- 


<0.10 


4.7 


SJ02 


(79983) 


1999 DF, 


3.3 


6.65/9.2 


0.40±0.02 


6.1 


LL06 


(40314) 


1999 KRi f , 


(5.840 or 5.929)±0.001 


(11.680 or 11.858)±0.002 


0.18±0.04 


5.8 


SJ02 


(47171) 


1999 TC 36 


6.21 


- 


0.06 


4.9 


O03b 






_ 


_ 


<0.07 




LL06 






- 


- 


<0.05 




SJ03 


(29981) 


1999 TD, 


7.71±0.02 


- 


0.65±0.05 


8.8 


O03b 






- 


15.45 


0.65 




R03.C03 








15.382 


0.41±0.08 




M04 






- 


15.3833 


0.53±0.03 




R05a 






5.8 


- 


0.65±0.05 




COO 




1999 TZ, 


_ 


10.438 


<0.1 


9.3 


M08 






5 


- 


0.15±0.02 




D08 






5.211 


10.422 


0.06±0.01 




Th09 


(80806) 


2000 CM 105 


- 


- 


<0.14 


6.3 


LL06 




2000 CP I04 


- 


- 


0.06 


6.7 


R08 




2000 FV53 


3.79 


- 


0.07 ±0.02 


8.2 


TB05 






- 


7.5 


0.07 




TB05 




2000 CP104 


- 


- 


0.06 


6.7 


R08 


(47932) 


2000 GN171 


- 


8.329±0.05 


0.61±0.03 


6.0 


SJ02 






- 


8.329±0.05 


0.60±0.03 




D08 


(82075) 


2000 YW, 34 


- 


- 


<0.10 


5.0 


SJ03 


(150642) 


2001 CZ 31 




4.71/5.23 


0.21±0.02 


5.7 


LL06 






- 


- 


<0.20 




SJ02 


(82158) 


2001 FP 185 


- 


- 


<0.06 


6.1 


SJ03 


(82155) 


2001 FZ m 


- 


- 


<0.06 


6.2 


SJ02 




2001 KD77 


- 


- 


<0.07 


5.8 


SJ03 




2001 QF 2 98 


- 


- 


<0.12 


4.7 


SJ03 


(139775) 


2001 QG298 


6.89±0.0002 


13.7744±0.0004 


1.14±0.04 


7.0 


SJ04 


(8861 IB) 


2001 QT 2 , 7B 


4.75 


- 


0.6 


5.5 


Os03 






5.50±0.02 


- 


- 




K06 


(42301) 


2001 UR163 


- 


- 


<0.08 


4.2 


SJ03 


(126154) 


2001 YH140 


6.22/8.45±0.05/12.99 


12.38/16.80/26.72 


0.16±0.05 


5.4 


O06 






- 


13.25 


0.21±0.04 




S07 






6.19/8.40/13.36 


12.38/16.80/26.72 


0.16±0.05 




Th09 


(55565) 


2002 AW„ 7 


6.49/8.87±0.01/13.94/15.82 


17.74 


0.04±0.01 


3.3 


O06 






8.78 


17.56 


0.04±0.01 




Th09 


(95626) 


2002 GP 32 


-3.3 


-6.6 


> 1 


6.7 


K06 






- 


- 


<0.03 




S07 


(95626) 


2002 GZ 32 


- 


5.80±0.03 


0.15±0.03 


6.8 


D08 


(73480) 


2002 PN 34 


(4.23 or5.11)±0.03 


(8.45 or 10.22)±0.06 


0.18±0.04 


8.2 


O03b 


(55636) 


2002TX 300 


(8.12 or 12.10)±0.08 


(16.24 or 24.20)±0.08 


0.08±0.02 


3.3 


SJ03 






7.89±0.03 


15.78 


0.09±0.08 




O04 






4.08 


8.16 


0.04±0.01 




Th09 


(55637) 


2002 UX 25 


7. 19 or 8.4 


14.38 


0.21 ±0.06 


3.6 


R05b 






- 


16.782 


0.13±0.09 




R05a 






- 


- 


<0.06 




SJ03 


(555638) 


2002 VE95 


6.76/6.88/7.36/9.47 


- 


0.08±0.04 


5.3 


O06 






- 


- 


<0.06 




SJ03 






4.90/6.18/9.97/17.23 


9.80/12.36/19.94/34.46 


0.04±0.01 




Th09 


(136204) 


2003 WL 7 


8.24 


16.48 


0.04±0.01 


8.7 


Th09 




2005 CB79 


6.76 


13.52 


0.13±0.02 


5.0 


Th09 


(145451) 


2005 RM 43 


6.71 


13.42 


0.04±0.01 


4.4 


Th09 


(145452) 


2005 RN 43 


5.62 or 7.32 


11.24 or 14.64 


0.04±0.01 


3.9 


Th09 


(145453) 


2005 RR 43 


7.87 


15.74 


0.06±0.01 


4.0 


Th09 


(145486) 


2005 UJ 438 


4.16 


8.32 


0.11±0.01 


10.5 


Th09 




2007 UL, 26 or 2002 KY 14 


3.56 or 4.20 


7.12 or 8.40 


0.10±0.01 


9.4 


Th09 


(120347) 


2004 SB 60 


6.09 or 8.01 


12.18 or 16.02 


0.03±0.01 


4.4 


Th09 


(144897) 


2004 UX 10 


5.68 


11.36 


0.08±0.01 


4.7 


Th09 


(119979) 


2002 WC19 






<0.03 


5.1 


S07 




2003 AZ 84 


4.32/5.28/6.72/6.76±0.01 




0.14±0.03 


3.6 


O06 








5.72±0.05 


0.14±0.03 




SJ02 






6.79 


13.58 


0.10±0.01 




Th09 


(120061) 


2003 CO, 


4.51 


9.02 


0.06±0.01 


8.9 


Th09 




2003 BF 91 


9.1/7.3 




1.09±0.25 


11.7 


TB05 




2003 BG91 


4.2/4.5/4.6/4.9 




0.18±0.75 


10.7 


TB05 




2003 BH91 






0.42 


11.9 


TB05 






2.8 




<0.15 




TB05 
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Object 


Designation 


Single Peak 
Period [h] 


Double Peak 
Period [h] 


Amplitude 
[mag.] 


Absolute 
magnitude' 1 


Ref b 


(120132) 


2003 FY 128 


- 


- 


<0.08 


5.0 


S07 






>7 


- 


- 




DOS 






8.54 


17.08 


0.12±0.01 


5.0 


Th09 


(120178) 


2003 OP32 


4.845±0.03 


- 


0.26±0.04 


4.1 


Ra08 






4.05 


8.10 


0.10±0.01 




Th09 




2003 QB,, 2 


- 


- 


0.46 


7.0 


R08 




2003 QY 9(M 


3.4±1.1 




0.34±0.06 


6.3 


K06a 




2003 QY90B 


7.1±2.9 




0.90±0.18 


6.3 


K06a 




2003 QYni 






0.60 


6.6 


R08 


(174567) 


2003 MW 12 


5.90/7.87 


11.80/15.74 


0.06±0.01 


3.6 


Th09 


(84922) 


2003 VS 2 


3.71 or 4.39 


7.42 


0.17±0.01 


4.2 


O06,Th09 








7.41±0.01 


0.21±0.02 




S07 


(90568) 


2004 GV 9 






<0.1 


4.0 


S07 








5.86±0.03 


0.16±0.03 




D08 


(120348) 


2004 TY 364 


5.85±0.01 


11.7±0.01 


0.22±0.02 


4.5 


S07 
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Table 2. Most important correlations found using the light curve 
parameters and orbital/physical variables as described in the text. 
The parameters presented are the light curve amplitude (Ampl.), 
the perihelion (q) and aphelion (Q) distances, inclination (i), ec- 
centricity (e) and absolute magnitude (H). 



Population 


Correlated Magnitudes 


P 


SL(%) 


n 


All KBOs 


Ampl. vs. Q 


-0.22 


96.96 


100 




Ampl. vs. H 


0.35 


99.94 


100 


Classical 


Period vs. B-V 


0.66 


98.49 


17 




Ampl. vs. Q 


-0.36 


97.95 


42 




Ampl. vs. H 


0.49 


99.80 


42 




Ampl. vs. i 


-0.37 


98.16 


42 




Ampl. vs. e 


-0.34 


96.74 


42 


Hot 


Period vs. B-V 


0.66 


98.11 


15 


Cold 


Period vs. Q 


-0.80 


94.97 


7 


Centaurs 


Period vs. B-V 


0.49 


91.40 


14 




Period vs. e 


0.47 


94.04 


18 


Plutinos 


Period vs. H 


-0.48 


94.14 


17 


SDO's 


Period vs. V-R 


-0.65 


91.55 


8 




Ampl. vs. H 


0.57 


92.86 


14 
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Table 3. Some results using the simple model. Listed are the minimum and maximum rotational period that a Jacobi ellipsoid can 
rotate. The percentage of different kinds of objects (Jacobi/MacLaurin/non-equilibrium) are presented. 



Density P [min] P [max] AP P [eq] % Jacobi % MacLaurin % No Eq. 



kg/m 3 


hrs 


hrs 


hrs 


hrs 










Q QQ 

7.77 


If) Q1 


n Q9 


8 88 
o.oo 


U.J 1 


99 f)4 


/ 1 .^rj 




8 88 
o.oo 


Q QQ 

7.77 


1 1 1 

1.11 


7 QQ 

/ .77 


Q ZL8 

7.^fO 


9fi f\l 


QPl 


600 


8.27 


8.88 


0.61 


7.27 


7.53 


36.34 


56.12 


700 


7.50 


8.27 


0.77 


6.66 


10.33 


41.29 


48.37 


800 


7.05 


7.74 


0.69 


6.31 


10.38 


46.26 


43.35 


900 


6.66 


7.27 


0.61 


5.85 


10.28 


53.52 


36.19 


1000 


6.31 


7.05 


0.74 


5.58 


12.61 


55.63 


31.75 


1200 


5.71 


6.48 


0.77 


5.10 


14.09 


62.11 


23.79 


1500 


5.10 


5.71 


0.61 


4.61 


11.92 


72.31 


15.76 


2000 


4.44 


5.00 


0.56 


4.00 


10.04 


82.77 


7.19 


2500 


4.00 


4.44 


0.44 


3.58 


6.72 


90.19 


3.08 


3000 


3.63 


4.06 


0.43 


3.24 


4.92 


94.06 


1.02 



